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The grand canonical potential for a static quark-antiquarkpair at non- vanishing chemical potential 
was determined at several temperatures. The dynamical staggered simulations were carried out 
with 2+1 physical quarks along the Lines of Constant Physiscs (LCP). 
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1. Introduction 

Forces among infinitely heavy quarks separated by distance r are of significance both at T = 
and when they are surrounded by an interacting medium. This paper is an extension of [1], 
where a reweighting technique was used to compute the grand canonical potential of a qq pair 
suiTounded by a medium of finite chemical potential /i along the transition line. In this work two 
other temperatures, one above and one below the transition temperature, will be examined. All 
simulations were carried out for 2+ 1 flavours with physical quark masses along a line of constant 
physics (LCP). 

2. Line of constant physics 

As we want to perform all our finite temperature simulations on lattices of size 4 x 12 3 , we 
have to control the temperature T = 1/ (4a) by the lattice spacing a. But on changing a it is also 
necessary to tune the mass parameters m u j and m s accordingly, where the hat denotes the quantities 
in lattice units. To determine this so called line of constant physics, we carried out zero temperature 
simulations for several values of the gauge coupling /3 and of the mass parameters (see Table 1). 
We used 2 + 1 flavours with dynamical staggered fermions. A few hundred configurations were 
generated for each parameter set using the R-algorithm and setting the microcanonical stepsize to 
half of the light quark mass. We measured m K , m p , niK, fo and a at 3 quark masses for each /3 in 





m u ,d 


m s 


5.09 


0.020, 0.040, 0.060 


0.250 


5.16 


0.020, 0.040, 0.060 


0.250 


5.19 


0.020, 0.040, 0.060 


0.250 


5.32 


0.010, 0.020, 0.030 


0.168 


5.42 


0.010, 0.020, 0.030 


0.126 


5.50 


0.008, 0.016, 0.024 


0.100 


5.57 


0.008, 0.016, 0.024 


0.100 



Table 1: Parameters used for the zero temperature simulations. 

order to extrapolate to physical masses. The lattice sizes ranged from 24 x 12 3 for coarse lattices 
over 32 x 16 3 to 48 x 24 3 for finer lattices, so that even for the lightest quark masses m n L > 4, 
where L is the spatial extent of the lattice. The ratio m-^/mp was set to its physical value (0. 179) in 
order to fix m u j as a function of j3 and m n /niK was set to its physical value (0.278) to determine 
rhu,d/>n s . The values of m uc i and m s as the function of j3 are shown in Figure [j]. A weighted average 
of m p , ro and a was used to set the scale (Figure ||). We extrapolated the ratio of the strange and 
light quark masses to the continuum limit (assuming 0(a 2 ) scaling) and got 

fit 

—5- =24.2(7) , (2.1) 

which is in agreement with chiral perturbation theory [2]. A similar line of constant physics using 
3 flavours is presented in [3]. 
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Figure 1: 20 m u j and m s as a function of /3 so 
that the m K /m p and m K /m.K are kept constant. 
The black points denote the parameters used 
in the finite temperature simulations. 

3. Reweighting technique 

As the fermion determinant is no longer positive for non-zero chemical potential any kind 
of importance sampling is hindered. We used the overlap-improving multiparameter reweighting 
technique [4]. The expectation value of an operator & at non- vanishing chemical potential can be 
rewritten as the ratio of two expectation values at }Xq = 0, where importance sampling can be used 
to generate the configurations. 

JdW 0(0 JdW &(o' [detMGu )]" //4 e~ Ss(ft) 



Figure 2: The scale I /a as a function of j3. 
The black points denote the parameters used 
in the finite temperature simulations. 



^ fd&(0 r^w'[detMGu )] n ' /4 e- s *(ft) (^'W 

[detM(n)] n f/ 4 e- s *W 



(3.1) 



where ft) = [detM(ju)] 



n//4 e -VP). 



ft) 



[detM^o)]"' 74 ^^ ' 

The details of the determinant evaluation can be found in [5]. As shown by eq. ( |3~l| ) a reweight- 
ing in other parameters like j8 is also possible and it has the advantage that it can maximize the 
overlap or trace out the phase transition line. There are however two limitations to the reweighting 
approach. First the overlap between the sample generated at /x = and the target sample at jX ^ 
decreases with growing fx. One can define a quantity which describes the overlap in several ways 
[6]. In this work it is defined as the ratio of the number of configurations in two different samples 
#(set')/#(set). Here (set) marks the sample generated at }X = and (set') the sample obtained after 
an accept-reject step with probability min(l, | Ct>/ /ct>/_ 1 1 ) - The index i stands for the ith configura- 
tion in (set). This overlap can then be maximized by searching for the optimal path in the jtx-jS 
plane, as shown in Figure ||. The line of maximal overlap for T around Tq matches the transition 
line, whereas for T = 0.9 Tq and T = 1.33 Tc it is approximately a line of constant temperature at 
least in the area where the reweighting is reliable. A second limitation, which is not measured by 
the overlap defined above is the sign problem. Its impact on the measurements is contained in the 
jackknife error. 
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Figure 3: Overlap in the plane around the transition line: The transition line (taken from [7]) matches 
the line of maximal overlap quite well, as expected [8]. 

4. Static qq free energy at T ^ and ji ^ 

We used the Polyakov-loop correlator to compute 0^(r,/i,r), which denotes the difference in 
the free energies with and without a static qq pair separated by distance r. 

e-^ T ^ r)/T+c(a) = {L{r)L\0)) T where L(r) = -tr ]J %(r,r) (4.1) 

Here C is a renormalization constant, which was fixed by first normalizing the zero temperature 
potential demanding V q q(ro) = and then removing this shift from the free energy at T ^ and 
jU 0, as it is the self-energy of the static quark pair at a given lattice spacing. Note that this 
type of renormalization is manifestly temperature independent. The finite temperature simulations 
were carried out using dynamical staggered quarks with n/ = 2 + 1 flavours on 4 x 12 3 lattices. 
The R-algorithm was used, the microcanonical stepsize was half of the light quark mass. We 
generated about 150000 configurations and calculated the determinant after every 50 trajectories 
giving us about 1500 independent configurations at T = 0.9T C and T = l.33T c and about 3000 at 
T = T c . Figure [|-^ show (j) q q for the three temperatures. The potential at = MeV is plotted 
in black (upper curves), that at /j.b = 270 MeV in red (lower curves). For all temperatures only 
slight changes can be observed between the potential at zero and finite jU. This is expected as the 
reweighting was performed along the lines of maximal overlap. For all temperatures the potential 
reaches an asymptotic value at large distances, due to the screening of the medium. Generally it 
can be stated, that for higher /I values the potential reaches a lower asymptotic value. Looking at 
the different scales of the axes one observes that the change in the potential increases with falling 
temperature, although being generally rather small. It can be observed that in the deconfined phase 
the potential flattens out at shorter distances. This is also demonstrated in Figure ||, where in blue 
(lowest curve) the potential at Hb = 270 MeV reweighted not along the transition line but along the 
line of constant temperature is drawn. Here the reweighting line cuts into the deconfined phase. 
In Figure |7| the free energy at a fixed distance is plotted as a function of the baryonic chemical 
potential. Particulary in the confined phase and also along the transition line one observes an 
increase of the errorbars due to the sign problem. A similar study using a Taylor expansion in 
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(jl/T) with 2 degenerate quark flavours and masses corresponding to m % « 770 MeV [11] was 
also presented at the symposium. A thorough comparison is difficult, but we observe that the shifts 
[(j) q q(T,0,°°) — <j>qq(T,n,°°)] are in general agreement, although our values tend to be about 15% 
smaller. 

5. Conclusions 

We have calculated the grand canonical potential of a static qq pair for three different lines on 
the T-fL plane: along the transition line as well as along one line above and along one line below 
it. All computations were done at physical quark masses on the line of constant physics using 
an overlap improving multiparameter reweighting in the gauge coupling /3 and in the chemical 
potential }i. We observed only a small change in the potential with growing /I in the parameter 
range where our analysis was earned out. 

Acknowledgments The computations were carried out at Edtvos University on the 330 processor PC 
cluster of the Institute for Theoretical Physics and the 1024 processor PC cluster of Wuppertal University, 
using a modified version of the publicly available MILC code [9] and a next-neighbour communication 
architecture [10]. 
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Figure 4: The grand canonical potential for 
T = 0.9 Tq'. The upper curve (black) denotes 
the potential at fj,g = MeV, whereas the 
lower curve (red) is the potential reweighted 
along the best reweighting line to jig = 270 MeV. 
In order to reduce the noise and to make the plot 
more readable, a box averaging was used. The 
box size ranged from 0.008 fm, at short distances 
to 0.2 fm at large distances. 
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Figure 5: The grand canonical potential for T = 
1 .0 Tq\ In this case the upper curve (black) denoting 
the potential at fj,g = MeV and the slightly lower 
curve (red) denoting the potential reweighted along 
the best reweighting line lie nearly on top of each 
other, due to the reweighting along the best reweight- 
ing line. The lower curve (blue) is the potential 
reweighted along a line of constant temperature to 
Hb = 270 MeV. 
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Figure 6: The grand canonical potential 
for T = 1 .33 Tq. The upper curve (black) denotes 
the potential at jj,g — MeV, whereas the 
lower curve (red) is the potential reweighted 
along the best reweighting line to JJ-b = 270 MeV. 
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Figure 7: ^qq as a function of /ig for a fixed 
distance r = y/lO. 



178/6 



